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ABSTRACT 

We present an overview of four deep phase-constrained Chandra HETGS X-ray observations 
of S Ori A. Delta Ori A is actually a triple system which includes the nearest massive eclipsing 
spectroscopic binary, 6 Ori Aa, the only such object that can be observed with little phase- 
smearing with the Chandra gratings. Since the fainter star, 5 Ori Aa2, has a much lower X-ray 
luminosity than the brighter primary (i5 Ori Aal), 6 Ori Aa provides a unique system with which 
to test the spatial distribution of the X-ray emitting gas around 6 Ori Aal via occultation by 
the photosphere of, and wind cavity around, the X-ray dark secondary. Here we discuss the 
X-ray spectrum and X-ray line profiles for the combined observation, having an exposure time 
of nearly 500 ks and covering nearly the entire binary orbit. The companion papers discuss 
the X-ray variability seen in the Chandra spectra, present new space-based photometry and 
ground-based radial velocities obtained simultaneous with the X-ray data to better constrain the 
system parameters, and model the effects of X-rays on the optical and UV spectra. We find 
that the X-ray emission is dominated by embedded wind shock emission from star Aal, with 
little contribution from the tertiary star Ab or the shocked gas produced by the collision of the 
wind of Aal against the surface of Aa2. We find a similar temperature distribution to previous 
X-ray spectrum analyses. We also show that the line half-widths are about 0.3 — 0.5 times the 
terminal velocity of the wind of star Aal. We find a strong anti-correlation between line widths 
and the line excitation energy, which suggests that longer-wavelength, lower-temperature lines 
form farther out in the wind. Our analysis also indicates that the ratio of the intensities of the 
strong and weak lines of Fe XVII and Ne X are inconsistent with model predictions, which may 
be an effect of resonance scattering. 
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— stars: mass-loss — X-rays: stars 
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1. Introduction 

Massive 0-type stars, though rare, are a pri¬ 
mary drivers of the chemical, ionization, and pres¬ 
sure evolution of the interstellar medium. The 
evolution of these stars from the main sequence 
to supernova depends on their mass and is signifi¬ 
cantly affected by stellar wind mass-loss. Our best 
estimates of mass, radius, and luminosity for O 
stars come from direct dynamical analyses of pho¬ 
tometric and radial velocity variations in massive, 
eclipsing binaries. However, because massive stars 
are rare and massive binaries which have been 
studied in detail rarer still (of the 2386 systems 
listed in the Ninth Catalog of Spectroscopic Bina¬ 
ries, only 82 of them have 0 -type components), 
direct dynamical determinations of stellar param- 
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eters are only known for a few systems. 

Current uncertainties regarding the amount 
and distribution of mass lost through stellar winds 
are even larger, since it is difficult to deter¬ 
mine stellar wind parameters in a direct, model- 
independent way. Radiatively driven stellar winds 
have mass-loss rates of M ~ 10 “® — 10 “^ M 0 yr“^ 
(for a review, see Kudritzki & Puls 2000). How¬ 
ever, observationally determined mass-loss rates 
have been estimated, in many, if not most cases, 
using an idealized smooth, spherically symmetric 
wind. Stellar winds are probably not spherical; 
variations of photospheric temperature with lati¬ 
tude are inevitable because of stellar rotation (and 
tidal deformation of stars in binaries), and these 
temperature variations will produce latitudinally 


dependent wind densities and velocities (Owocki 


et al. 1996). Stellar winds are not smooth ei¬ 


ther; the radiative driving force is inherently un¬ 
stable to small velocity perturbations, and wind 
instabilities are expected to grow into dense struc¬ 
tures (clumps) distributed through the wind. In 
addition, clumps can also be produced by sub¬ 
surface convective zones in massive stars caused 
by opacity peaks associated with the ionization 
state of helium and iron (Cantiello et al. 2009). 
Wind clumps play an important role in determin¬ 
ing the overall mass-loss rate, since they carry 
most of the mass but occupy little volume. An 
outstanding question is to determine the number 
and mass/spatial distribution of embedded wind 
clumps. 

Collisions between clumps, or between clumps 
and ambient wind material at high differential ve¬ 
locities can produce pockets of hot shocked gas 
embedded in the wind. Given wind speeds of up 
to thousands of kilometers per second, these em¬ 
bedded wind shocks should generate observable 


X-ray emission (as originally proposed by Lucy 


& White 1980). There have been efforts to de¬ 


termine the fraction of the wind that is clumped, 
and the radial distribution of the embedded wind 
shocks, through analysis of the X-ray radiation 
they produce. High spectral resolution X-ray grat¬ 
ing spectrometry provides a unique tool to de¬ 
termine the properties of the X-ray emitting hot 
shocked gas produced by embedded wind clumps. 
In particular, the forbidden-to-intercombination 
{f /i) line ratios of strong He-like transitions, and 
analysis of profiles of H-like ions and other strong 
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lines from high resolution spectra (mostly from the 
Chandra and XMM grating spectrometers) indi¬ 
cate that significant X-ray emission exists within 


1 to 2 radii of the stellar photosphere (Waldron 

& Cassinelli 

2001 

Leutenegger et al. 2006[ jWal- 

dron & Cassinelli|2007 

). X-ray lines of strong Lya 


transitions (mainly O VIII, Ne X, Mg XII, Si XIV, 
and S XVI) show prohles ranging from broad and 
asymmetric to narrow and symmetric, apparently 


dependent on stellar spectral type (Walborn et al. 


20091. Observed line profile shapes are an impor¬ 


tant probe of the radius of the maximum X-ray 
emissivity, modified by absorption from the over- 
lying, cooler, clumped wind. 

dumping-corrected mass-loss rates derived 
from the analysis of resolved X-ray emission lines 
( Oskinova et al.|2006 1 are generally in good agree¬ 
ment with predictions of line-driven wind theory, 
while mass-loss rates derived from analyses of re¬ 
solved X-ray emission lines are lower (by a factor 
of a few) if clumping is not taken into account 
(Cohen et al. 2014). Reducing mass-loss rates 


by such a large factor would significantly influ¬ 
ence our understanding of the ultimate evolution 
of massive stars. However, while important wind 
properties, such as the onset radius of clumping, 
the fraction of the wind that is clumped, and the 
radial distribution of clumps through the wind, 
have been indirectly inferred from detailed X-ray 
line analysis ( Oskinova et al.|[2006 Owocki & Co¬ 
hen 2006 Herve et al. 2013), to date, there have 


been no attempts to determine these properties 
directly. In this paper, we try to directly con¬ 
strain the location of the X-ray emitting gas in 
the wind of a massive eclipsing binary, 5 Ori Aa, 
via occultation by the companion star of the hot 
gas embedded in the primary’s wind. 

Delta Ori (Mintaka, HD 36486, 34 Ori) is a vi¬ 
sual triple system composed of components A, B, 
and C. Delta Ori A itself is composed of a massive, 
short period close eclipsing system 6 Ori Aa, and 
a more distant component, S Ori Ab, which orbits 


6 Ori Aa with a period of 346 years (Tokovinin 


et al. 2014). The inner binary, 6 Ori Aa, is the 
nearest massive eclipsing system in the sky. It 
consists of a massive 09.5 H primary (star Aal) 
-I- a fainter secondary (star Aa2, B2V-B0.5 HI), 
in a high-inclination {i > 67°), short period jP = 
5°*.7324), low eccentricity (e « 0.1) orbit (Hart- 


|1904[[Stebbins|19T5l|Koch &: Hriviiak 1981 


mann 


Harvin et al. 2002 Mayer et al. 2010). Because 


it is nearby, bright, with a high orbital inclina¬ 
tion, S Ori Aa is an important system since it 
can serve as a fundamental calibrator of the mass- 
radius-luminosity relation in the upper HR dia¬ 
gram. It is disconcerting, though, that published 
stellar masses for the primary star 6 Ori Aal are 


different by about a factor of two (Harvin et al. 
2002| [Mayer et al.||2010r 


Delta Ori Aa is also a bright X-ray source (Long 


& White 1980 Snow et al.j |1981| jCassinelli &: 
Swank|19^ ) and is the only eclipsing short-period 


0-type binary system that is bright enough to be 
observable with the Chandra gratings with little 
phase smearing, offering the chance to study of 
variations of the X-ray emission line profiles as a 
function of the orbital phase. 


Since the luminosity of the secondary, 5 Ori Aa2, 
is less than 10% that of the primary, and since 
X-ray luminosity scales with stellar bolometric lu- 


minosity (Pallavicini et al.|1981 

Chlebowski et al. 

1989 

Berghoefer et al.|1997 

) for stars in this mass 


range, it should also be less than 10% as bright in 
X-rays as the primary. Thus the X-ray emission 
from the system is dominated by the hot gas in 
the wind of the primary star. Therefore, occulta¬ 
tion of different X-ray-emitting regions in the wind 
of 6 Ori Aal by the photosphere and/or wind of 
the X-ray faint secondary, 6 Ori Aa2, presents 
the opportunity to directly study the radial dis¬ 
tribution of the hot shocked gas in the primary’s 
wind, by measuring occultation effects in X-ray 
line emission as a function of ionization potential 
and orbital phase. Since X-ray lines of different 
ionization potentials are believed to form at differ¬ 
ent radial distances above the primary’s surface, 
differential variations in the observed set of X- 
ray lines as a function of orbital phase allow us 
to probe the hot gas distribution within the pri¬ 
mary wind’s acceleration zone, where most of the 
X-ray emission is believed to originate. He-like 
ions in the X-ray spectrum provide a complemen¬ 
tary measure of the radial distribution of the hot 
gas, since these lines are sensitive to wind den¬ 
sity and the dilute ambient UV field. This makes 
5 Ori Aa a unique system with which to constrain 


^ Some progress has been recently made by Harmanec et al. 
l |2013^ and by [Richardson et al.| l |2015| l in disentangling 
lines of 5 Ori Aa2 from S Ori Aal and 5 Ori Ab in the 
composite spectrum 
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directly the spatial distribution of X-ray emitting 
clumps embedded in the wind of an important O 
star. The main challenge, however, is the rela¬ 
tively small size of S Ori Aa2 compared to the 
size of the X-ray emitting region, since the hot 
gas is expected to be distributed in a large volume 
throughout the stellar wind. 

This paper provides an overview of the X- 
ray grating spectra obtained during a 479 ksec 
Chandra campaign on 6 Ori Aa-|-Ab in 2012. 
The purpose of this project was to obtain high 
signal-to-noise observations with the Chandra 
High Energy Transmission Grating Spectrometer 


2. Stellar And System Parameters 


(HETGS; Canizares et al. 20051 of 6 Ori Aa over 


almost an entire binary orbit, including key orbital 
phases, with coordinated ground-based radial ve¬ 
locity monitoring at Ha and He I 6678 (primar¬ 
ily obtained by a group of amateur astronomers), 
and high precision, simultaneous photometry from 
space by the Canadian Space Agency’s Microvari¬ 
ability and Oscillations of Stars telescope {MOST, 
Walker et al. 2003). This paper provides an 


overview of the combined HETGS spectrum from 
our four observations, and is organized as follows. 
In Section we present a summary of the four 
observations and discuss the acquisition and re¬ 
duction of the data sets. Section presents an 
analysis of the zeroth-order image of the system 
to constrain the X-ray contribution of S Ori Ab to 
the observed X-ray emission. Section presents 
the temperature distribution and overall proper¬ 
ties of the strong emission lines in the combined 
spectrum of the four observations. Section dis¬ 
cusses the possible influence of the collision of the 
wind from the primary with the weak wind or pho¬ 
tosphere of the secondary, and the influence of any 
such collision on the wind’s thermal and density 
structure. We present conclusions in Section 
A series of companion papers presents the results 
of the variability analysis of the X-ray continuum 
and line emission (Nichols et al. 2015, in press. 


Paper H), the ground-based radial velocity and 
MOST space-based photometric monitoring and 


analysis (Pablo et al. 2015 in press. Paper HI), 
and a complete non-LTE analysis of the spectral 
energy distribution of 5 Ori Aa-l-b from optical 
through X-rays (Shenar et al. 2015 in press, Pa¬ 
per IV). 


The stellar parameters given by |Harvin et al.| 
(2002) and Mayer et al. (2010) differ significantly. 


and this difference has important consequences for 
our understanding of the evolutionary state of the 
system, and the influence of mass-loss and/or non¬ 
conservative mass transfer. Harvin et al. (2002) 


derived masses of M^ai = 11 . 2 M 0 and MAa 2 = 
5.6Mq for the primary and secondary stars, mak¬ 
ing the primary significantly overluminous for its 
mass (or undermassive for its spectral type). The 
radial velocity and photometric analysis of [Mayer] 
et al. ( 2010 ) were consistent with a substantially 
higher mass for the primary, MAai = 25Mq, af¬ 
ter a correction for perceived contamination of 
the radial velocity curve by lines from 6 Ori Ab. 
Whether the 09.5 H primary has a normal mass 
and radius for its spectral type is important for 
understanding the history of mass exchange/mass- 
loss from S Ori Aa, and how this history is related 
to the current state of the radiatively driven wind 
from the primary. 

An important goal of our campaign is to de¬ 
rive dehnitive stellar and system parameters for 
S Ori Aa. To this end, we obtained high-precision 
photometry of the star with the MOST satellite, 
along with coordinated ground-based optical spec¬ 
tra to allow us to obtain contemporaneous light- 
and radial-velocity curve solutions, and to dis¬ 
entangle the contributions from Aa 2 and/or Ab 
from the stellar spectrum. We also performed 
an analysis of the optical and archival lUE UV 
spectra using the non-LTE Potsdam Wolf-Rayet 


(PoWR) code (Grafener et al. 2002 Hamann & 
Grafener 2003). The light curve and radial ve¬ 


locity curve analysis is presented in Pablo et al. 
(2015), while the non-LTE spectral analysis is pre¬ 
sented in Shenar et al. (2015). Table summa¬ 
rizes these results. In this table, the values and 
errors on the parameters derived from the MOST 
photometry and radial velocities are given for the 
low-mass solution provided in Pablo et al. (2015). 
Note that we find better agreement between the 
derived stellar parameters (luminosities, masses, 
radii, and temperatures) and the spectral type 
of 6 Ori Aal if we use the cr-Orionis cluster dis¬ 


tance {d = 380 pc, Caballero & Solano 2008) for 


<5 Ori A, rather than the smaller Hipparcos dis¬ 
tance. Therefore, we adopt D = 380 pc as the 


4 


















































distance to S Ori A (for a full discussion of the 
distance to 6 Ori A, see Shenar et al. 2015). The 
spectral type of S Ori Aa2 is not well constrained; 


Harvin et al. (2002) assign it a spectral type of 


BO.5 III, while Mayer et al. (2010) do not assign 


a spectral type due to the difficulty in identifying 
lines from the star. Shenar et al. (2015) assign an 
early-B dwarf spectral type to 6 Ori Aa2 (« BIV). 

3. Previous X-ray Observations 

X-ray emission from S Ori was first tentatively 


identified via sounding rocket observations (Fisher 


& Meyerott 1964). X-ray imaging spectrometry 
of 6 Ori A at low or modest resolution was ob¬ 


tained by the EINSTEIN (Long & White 1980), 
ROSAT dHaberl fc White|1993D , and ASCA ( |Cor- 
coran et al. 1994) X-ray observatories. Its X-ray 
luminosity is typically ~ gj.gg g-i^ with 

Lx|L^|nl « 10“^ in accord with the canonical re- 


lation for massive 

stars 

(Pallavicini et al. 

1981 

Ghlebowski et al. 

1989 

Berghoefer et al. 

1997). 


The X-ray spectrum of 5 Ori A was observed at 
high resolution by X-ray grating spectrometers on 
Chandra in two previous observations at restricted 
orbital phases. An analysis of a fifty kilosecond 
Chandra HETGS spectrum from 2000 January 13 


by Miller et al. (2002) revealed strong line emission 


from O, Ne, Mg, and Fe, along with weaker emis¬ 
sion from higher-ionization lines like Si XIII and 
S XV, and unusually narrow line half-widths of 
« 400 km s“^. Using a simple analysis taking into 
account dilution of the photos pheric UV field and 
a 1 /r^ falloff in wind density, [Miller et al. (2002) 
derived formation regions for the dominant He-like 
ions Mg XI, Ne IX, and O VII extending just above 
the stellar photosphere to 3-10 times the photo- 
spheric radius. An analysis of a 100 ks Chandra 
Low Energy Transmission Grating Spectrometer 
(LETGS; Brinkman et al. 1987) -I- High Resolu¬ 
tion Camera observation from 2007 November 09 


by Raassen & Pollock (2013) also showed that the 


Mg XI, Ne IX, and O VII emission regions ex¬ 
tend from 2-10 stellar radius, and showed that the 
longer wavelength ions like N VI and C V form at 
substantially greater distances from the star (50- 
75 times the stellar radius), and that the spectrum 
could be modeled by a three-temperature plasma 
in collisional ionization equilibrium with temper¬ 
atures of 0.1, 0.2, and 0.6 keV. 


4. New Chandra Observations 

A listing of the Chandra observations of S Ori Aa-|-Ab 
obtained as part of this campaign is given in Ta¬ 
ble H These observations were obtained with the 
Chandra HETGS-I-ACIS-S spectrometric array. 

The HETGS consists of 2 sets of gratings: the 
Medium Resolution Grating (MEG) covering the 
range 2.5-26A and and the High Resolution Grat¬ 
ing (HEG) covering the range 1.2-15A; the HEG 
and MEG have resolving powers of A/AA « 1000 
at long wavelengths, falling to ^ 100 near 1.5A 


(Canizares et al. 2005). Eour observations cover¬ 


ing most of the orbit were obtained within a 9-day 
timespan to reduce any influence of orbit-to-orbit 
X-ray variations, for a combined exposure time 
of 479 ks. Table lists the start and stop HJD, 
phases, and exposure durations for the four in¬ 
dividual observations. Figure shows the time 
intervals of each observation superposed on the 
simultaneous MOST optical light curve of 6 Ori A 
(Pablo et al. 2015). The Chandra observations 


provide both MEG and HEG dispersed first order 
spectra as well as the zeroth order image. Due to 
spacecraft power considerations as well as back¬ 
ground count rate issues, it was necessary to use 
only five ACIS CCD chips instead of six; thus, chip 
S5 was not used. This means that wavelengths 
longer than about 19 A in the MEG plus-side 
dispersed spectrum and about 9.5 A in the HEG 
plus-side dispersed spectrum are not available. 
Therefore, the strong O VII line at 21 A was only 
observed in the MEG-1 order. The buildup of con¬ 
taminants on the ACIS-S optical blocking filters 
with time further degraded the long wavelength 
sensitivity for all first-order spectra. Each of the 
four observations experienced a large variation 
in focal plane temperature during the observa¬ 
tion. While a temperature-dependent calibration 
is applied to each observation in standard data 
processing, the calibration is based on a single 
temperature measurement taken at the end of the 
observation. In particular, the focal plane temper¬ 
ature for portions of each observation exceeded the 
temperature at which the temperature-dependent 
effects of charge transfer inefficiency (CTI) are 
calibrated (Grant et al. 2006). This could cause 


residual errors in the correction of pulse heights 
for those portions of the observations in the high- 
temperature regime. 
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Table 1: Stellar, Wind and System parameters for 5 Ori Aal+Aa2 from Analysis of the Optical, UV and 
X-ray spectra (Shenar et al. 20151 and the Solution to the MOST Light Curve and Ground-Based Radial 
Velocities ( Pablo et al.||2015 ). 


Parameter 

Method 

POWR Analysis® 

light curve & RV solution*’ 

Teff[kK] (Aal) 

Teff[kK] (Aa2) 
i?[R 0 ] (Aal) 

R[Rq] (Aa2) 

M[Mq] (Aal) 

M[Mq] (Aa2) 

L [log Lq] (Aal) 

L [log Lq] (Aa2) 

29.5 ±0.5 

25.6 ±3 

16.5 ± 1 

6-51L 

24:1:1° 

8.4® 

5.28 ±0.05 

4.2 ±0.2 

30 (adopted) 

24.11°;^7 

15.1 

5.0 

23.8 

8.5 

5.20 

3.85 

Uoo [km s“^] (Aal) 

2000 ± 100 


Voo [km s“^] (Aa2) 

1200 ® 


logM[M 0 /yr] (Aal) 

-6.4 ±0.15 


logM[M 0 /yr] (Aa2) 

< - 6.8 


Eb-v (ISM) 

0.065 ± 0.002 


Av (ISM) 

0.201 ±0.006 


logiVn (ISM) 

20.65 ±0.05 


Pld] 


5.732436®* 

Eq (primary min, HJD) 


2456277.790 ±0.024 

To (periastron, HJD) 


2456295.674 ±0.062 

a[i? 0 ] 


43.1 ± 1.7 

i [deg.] 


76.5 ±0.2 

uj [deg.] 


141.3 ±0.2 

uj [deg. yr“^] 


1.45 ±0.04 

e 


0.1133 ± 0.0003 

7 [km s“^] 


15.5 ±0.7 

Sp. Type (Aal) 

09.511“’®’®* 

Sp. Type (Aa2) 

B1V“ 

D [pc] 

380 (adopted) 


Notes: 


^Shenar et al. (2015); ^from the low-mass model solution of Pablo et al. (2015); 


Mai'z Apellaniz et al. 


I 20131; '^Mayer et al. (2010); ^Adopted assuming a spectral type of BIV. 
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Table 2: New Chandra Observations of 6 Ori Aa+Ab 


ObsID 

Start 

HJD 

Start 

Phase 

End 

HJD 

End 

Phase 

Midpoint 

HJD 

Midpoint 

Phase 

AT 

Days 

Exposure 

s 

Roll 

deg. 

14567 

2456281.21 

396.604 

2456282.58 

396.843 

2456281.90 

396.724 

1.37 

114982 

345.2 

14569 

2456283.76 

397.049 

2456285.18 

397.297 

2456284.47 

397.173 

1.42 

119274 

343.2 

14570 

2456286.06 

397.450 

2456287.52 

397.705 

2456286.79 

397.578 

1.46 

122483 

83.0 

14568 

2456288.67 

397.905 

2456290.12 

398.159 

2456289.39 

398.032 

1.45 

121988 

332.7 


Each ObsID was processed using the stan¬ 
dard processing pipeline used in production of 
the Chandra Transmission Grating Data Archive 
and Catalog (TGCAT; |Huenemoerder et al. 20111. 
Briefly, event filtering, event transformation, spec¬ 
tral extraction, and response generation are done 
with standard Chandra Interactive Analysis of Ob¬ 
servations software tools ( Fruscione et al.|20'06 ) as 
described in detail by Huenemoerder et al. (20111. 
This pipeline produces standard X-ray events, 
spectra, responses, effective areas, aspect his¬ 
tograms, and light curves. We used version 4.5.5 
of the Chandra Calibration Database (CALDB), 
along with CIAO version 4.5 & 4.6 in the analysis 
presented here. In order to examine variability, 
the data were also divided into ~ 10 ks segments, 
and spectra, response files, effective areas and light 
curves were generated for each segment. Analy¬ 
sis of the time-sliced data is presented in |Nichols| 


et al. (2015). 


5. Analysis of the X-ray Image 


The 6 Ori Aal,2 inner binary is orbited by a 
more distant tertiary component (5 Ori Ab) at 
a current projected separation of 0".3 with an or¬ 
bital period of « 346 years (Tokovinin et al. IMIf . 
This separation is just below the spatial resolu¬ 
tion of Chandra, and thus Chandra imaging ob¬ 
servations allow us to spatially examine the X-ray 
contribution from the Ab component. Figure 
shows unbinned zeroth-order images from our four 
HETGS-I-ACIS observations, along with the ex¬ 
pected location of Ab and the Aa pair at the times 
of the Chandra observations in 2012. 


To constrain the X-ray contribution of 5 Ori Ab, 
we generated zeroth-order images for the four in¬ 
dividual pointings listed in Table using the 
Energy-Dependent Subpixel Event Repositioning 
(EDSERp^ method to generate images with a 
pixel size of 0".125. We generated images in 0.3-1 
and 1-3 keV bands, but found no significant dif¬ 
ferences in any of the four observations when we 
compared the soft and hard band images. For each 
image, we then applied the CIAO tool SRCEX- 
TENT to calculate the size and associated uncer¬ 
tainty of the photon-count source image or using 
the Mexican Hat Optimization algorithrrp^ 

The results of the SRCEXTENT analysis are 
given in Table The derived major and minor 
axes of each image are equal and consistent with 
the Chandra point spread function, ^ 0.3". The 
peak of the image is consistent with the location 
of the Aa component, and is about a factor of two 
farther than the Ab component. We conclude that 
the peak positions of the zeroth-order images in- 


^^http: //cxc .harvard, edu/ciao4.4/why/acissubpix.htnil 
^■^http: //cxc. harvard.edu/ciao/ahelp/srcextent.html 
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Fig. 1.— Timings of the Chandra observations along with the MOST light curve. The images above the 
plot show the orientations of i5 Ori Aal and Aa2 near the midpoint of the observation according to the 
photometric and spectroscopic analysis of Pablo et al. (2015). In the images, the orbital angular momentum 
vector lies close to the plane of the paper and points to the top of the page. 


dicate that Aa is the primary X-ray source, with 
little or no contribution from Ab. Our analysis 
also suggests that the ObsID 14568 image may be 
slightly elongated, which may indicate a possible 
issue with the instrumental pointing or aspect re¬ 
construction for this observation. 

6. Combined Spectrum 

Figure shows the co-added spectrum from 
the four observations, with a total exposure of 
479 ks. This represents the second longest expo¬ 
sure yet obtained on a massive star at wavelengths 
< 8 A and a resolving power of A/AA > 400. The 
strongest lines are O VIII, Fe XVII, Ne IX & Ne X, 
Mg XI & Mg XII, and Si XIII. 


6.1. Temperature Distribution 


We modeled the combined spectrum with a 
combination of absorbed collisional ionization 
equilibrium models using the Interactive Spec¬ 


tral Interpretation System (ISIS; Houck & Deni- 


cola 2000). The model we applied includes two 


low-temperature components seen through a com¬ 
mon absorption component, plus a third hotter 
component with its own absorption component 
to account for any contribution from a hot col¬ 
liding wind region embedded within the wind of 
the binary (see Section below). In ISIS ter¬ 
minology, the mode we used was “(xaped(l) -I- 
xaped(2)) * TBabs(3) -I- xaped(4) * TBabs(5)”, 
where “xaped” represents emission from an op¬ 
tically thin plasma in collisional ionization equi¬ 
librium based on the ATOMDB atomic database 


version 2.0.2 (Smith & Brickhouse 2000 Foster 


et al. 2012), and “TBabs” represents interstellar 


absorption (Wilms et al. 2000a). Solar abundances 


were assumed for both the emission and absorp- 
















































Table 3: SRCEXTENT Analysis Results 



Band 

Major Axis 

Minor Axis 

PA 

Peak distance Aa 

Peak distance Ab 

ObsID 

keV 

arcsec 

arcsec 

deg. 

arcsec 

arcsec 

14567 

0.3-1 

0.34 

0.33 

83.3 

0.19 

0.40 


1-3 

0.32 

0.28 

83.8 

0.19 

0.42 

14569 

0.3-1 

0.32 

0.32 

32.1 

0.23 

0.44 


1-3 

0.29 

0.28 

27.6 

0.25 

0.47 

14570 

0.3-1 

0.32 

0.32 

136.9 

0.09 

0.35 


1-3 

0.26 

0.22 

48.3 

0.08 

0.34 

14568 

0.3-1 

0.51 

0.32 

35.9 

0.24 

0.41 


1-3 

0.48 

0.25 

31.2 

0.24 

0.42 



Fig. 2.— Unbinned images from the four ObsIDs listed in Table ObsIDs, left to right: 14567, 14568, 
14569, and 14570. The positions of Aa and Ab are shown by the full and dashed circles, respectively. 
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tion components!^ This model is an approxima¬ 
tion to the actual temperature distribution and 
absorption, but is the simplest one we found that 
adequately describes the observed grating spec¬ 
trum. We allowed for velocity broadening of the 
emission lines, with turbulent velocity broaden¬ 
ing constrained to be less than roughly twice the 
maximum wind terminal velocity, 3000 km s“^. 
We allowed the line centroid velocities of the three 
emission components to vary, but found that over¬ 
all the line centroids are unshifted in the combined 
spectrum. Figure [^compares the best-fit model to 
the data, while the model components are given 
in Table In this table, we also convert the 
derived turbulent velocities Vturb to equivalent 
line half-widths at half maximum, using O VIII, 
Ne X and Mg XII for the low-, medium-, and 
high-temperature components, respectively. 


The derived temperature distribution is similar 
to that found by Miller et al. (2002) in their study 
of the 2000 January HETGS spectrum, and by 


Raassen & Pollock (20131 in their analysis of an 
LETGS spectrum from 2007 November. In gen¬ 
eral, aside from the overall weakness of the forbid¬ 
den lines compared to the model spectrum (which 
assumes a low-density plasma with no UV pho¬ 
toexcitation), the overall distribution of emission 
line strengths, and the continuum, are described 
reasonably well by the model. We note, in reality, 
that this three-temperature model is a simplified 
representation of the actual emission measure dis¬ 
tribution with temperature. This multitempera¬ 
ture model mainly provides us with an adequate 
approximation of the local (pseudo-) continuum in 
order to improve line fitting and modeling. 


6.2. Emission Lines 


The observed X-ray emission lines in our 
S Ori A spectrum provide important diagnostic 
information about the phase-averaged state of the 
hot gas within the wind of the system, and, as 
we show below, this is dominated by the shocked 
gas embedded within the wind of S Ori Aal, with 
little contribution (if any) from gas heated by 
the shock produced by the collision of the wind 
from S Ori Aal with the wind or photosphere of 


Shenar et al. 


(2015 1 show that N and Si are slightly sub¬ 
solar, but these differences are not significant for our anal¬ 


ysis. 


S Ori Aa2. The analysis of the set of emission 
lines depends on choice of line profile, continuum 
level, and accounting for line blends. 

6.2.1. Gaussian Modeling 

To better account for blends and uncertainties 
in the continuum level, we performed a Gaussian 
fit to the strong lines, allowing flux, line width, 
and centroid velocity to vary. These fits, shown in 
Figure were done using the three-temperature 
fit given in Section |6.1| above to define the con¬ 
tinuum and amount of line blending. We set the 
abundance of the element to be measured to zero, 
with the abundances of other elements set to solar 
and other parameters (temperature, absorptions) 
fixed at the values given in Section This proce¬ 
dure is useful to account for line blends, in particu¬ 
lar, for the Ne X line at 12.132 A, which is blended 
with an Fe XVII line at 12.124 A. We assumed sim¬ 
ple Gaussian line profiles for the line to be fit, and 
fit for both the Lyal and Lya2 lines, with line 
widths and velocities fixed for both components, 
and the intensity ratio of the LyQ!2 to the Lyal 
line set to the emissivity ratio at the temperature 
of peak emissivity. We used the Cash statistic and 
ISIS to perform the fits, simultaneously fitting the 
HEG and MEG ±1 order spectrum from all four 
observations simultaneously. Table[^shows the re¬ 
sult of fits of the H-like Lya lines, plus the strong 
Fe XVII line at 15.014 A. In general, the Gaussian 
fits are poor (the reduced Cash statistic > 1.5) 
except for the weak Si XIV line, though the asym¬ 
metries in the bright lines are not very strong. All 
of the line centroids are near zero velocity, though 
the Ne X line is blue-shifted at about the 2-a level. 
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Table 4: Best-Fit to the Combined HETGS spectrum. 
(APED 1 -h APED 2 ) * 1 -t APED 3 


The adopted model is 


Component 

Parameter 

Value 


Ti (MK) 

1.25 

1 

EMi (10^5 cm-3) 

4.46 


(km s ) 

1313 


HWHM (km s-A 

1094 


T 2 (MK) 

3.33 

2 

EM 2 (10^5 cm-3) 

0.87 


Vturb ,2 (km s“A 

1143 


HWHM (km s-A 

953 

Absorption 1 

Nh,i (10^^ cm“^) 

0.14 


Ts (MK) 

9.11 

3 

EM 3 (10^5 cm-3) 

0.26 


Vturb ,3 (km s“A 

685 


HWHM (km s-i) 

574 

Absorption 2 

Nh,2 (10^^ cm“^) 

0.24 

fx (ergs cm“^ s“^) (observed, 1.7— 25 A) 


8.2 X 10-12 

Lx (ergs s“^) (observed, 1.7 — 25 A) 


1.4 X 10^2 

log Lx/L }jqi 


-6.73 

EM-weighted Average Temperature (MK) 


1.94 


Table 5: Gaussian Fits to the H-like lines, plus Fe XVII 



A 

Flux 

V 

HWHM 


A 

10“® ph. s”! cm-2 

km s”! 

km s”! 

0 VIII 

18.967 

219^?o 

_q+37 

^-33 

918+^^ 

Fe XVII 

15.014 

76;^ 


971/// 

NeX 

12.132 

lOtJ 

-102/11 

726/tl 

Mg XII 

8.419 

l+° 

-*^-0 

1 0+33 
-*-^-55 

M7/11 

Si XIV 

6.180 

A oc-l-0.05 

— 49+45 
^^-134 

5441116 
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Fig. 4.— The combined MEG+HEG spectrum of 6 Ori A (in black) with the 3-component fit (shown in 
red) given in Table The model spectra, which assume low density and do not include effects of UV 
photoexcitation, generally overestimate the strength of the forbidden lines and underestimate the strengths 
of the intercombination lines, especially at longer wavelengths, most notably at OVII. 
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Best Fit Model to OVIll 


Best Ht Model to FeXVIl 



Wavelength (Ang.) 



Best Fit Model to NeX 


Best Fit Model to MgXIl 




Wavelength (Ang.) 


Wavelength (Ang.) 


Best Fit Model to SIXIV 



Wavelength (Ang.) 


Fig. 5.— Top to bottom, left to right: O VIII; Fe XVII; Ne X; Mg XII; and Si XIV. The lines are plotted 
in the velocity range of —3000 km to +3000 km s“^. The best-fit Gaussian profile, and the continuum 
derived from the model parameters given in Table is shown in red. Note that while most of the Lya lines 
are adequately described by a symmetric Gaussian, the Fe XVII and Ne X lines are not as well fit by simple 
Gaussian profiles as the other lines. This may be due to the effects of non-uniform X-ray line opacity, as 
discussed in Section [6.2.21 
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We also measured the forbidden (z), intercom¬ 
bination {x + y), and resonance components (w) 
above continuum for each of the helium-like ions 
(O VII, Mg X, Ne IX, and Si XIII) by Gaussian 
fitting. As before, we used the three-temperature 
fit given in Section 6.1 above to define the local 
continuum near the line region. Although the in¬ 
dividual intercombination components {x + y) are 
unresolved in the HETGS spectra for all of the 
He-like ions, we included a Gaussian line for the 
X and y lines, but restricted the centroid velocity 
and line widths to be the same for both the x and 
y components. Because the forbidden, intercombi¬ 
nation and resonance lines can have different spa¬ 
tial distributions throughout the wind, we allowed 
the widths, centroids, and line fluxes of these lines 
to vary individually. The forbidden component 
of the O VII line is weak, and, in addition, this 
line was only observed in the MEG-1 spectrum 
arm because ACIS-S chip S5 was turned off due 
to spacecraft power constraints. To increase sig¬ 
nal to noise for the O VII forbidden line, and for 
the weak Si XIII and S XV triplets, we included 
data from the 2001 HETG and 2008 LETG obser¬ 
vations when fitting. Figurej^shows the fits to the 
He-like lines, and Table shows the results of this 
three-Gaussian component fitting, while Table 
shows the R = z/{x + y) and G = {x -\- y + z)/w 
ratios. 


Table 7: R and G ratios 


ion 

R = z/{x + y) 

G = {x + y z)/w 

0 VII 

0.04 ±0.01 

0.94 ±0.26 

Ne IX 

0.27 ±0.10 

1.44 ±0.65 

Mg XI 

0.96 ±0.36 

0.95 ±0.37 

Si XIII 

1.77 ±0.18 

0.90 ±0.12 

S XV 

3.88 ±2.86 

0.72 ±0.74 
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Best Ftl Model from OVILfit.sl 


Best Fit Model from NelX_fit.sl 




Best Fit Model from MgXIJit.sl Best Fit Model from SiXIlUit.sl 




Best Fit Model from SXV_frt.sl 



Wavelength (Ang.) 


Fig. 6.— Top to bottom, left to right: O VII; Ne IX; Mg XI; Si XIII; S XV. The best fit, using a model of 
4 Gaussian lines (w, a;, y, & z components) and the continuum derived from the model parameters given in 
Table 1^ is shown in red. 
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Table 6: Gaussian fits to the He-like lines 



Centroid Velocity (km 

s-i) 

HWHM (km s- 


ion 

w 

x + y 

z 

w 

x + y 

z 

OVII 

166 ± 19 

-194 ± 18 

-810 ±384 

761 ± 14 

826 ± 40 

160 ± 270 

Ne IX 

-146 ± 166 

-410 ±231 

441 ± 466 

849 ± 138 

1057 ±222 

1289 ± 49 

Mg XI 

8 ±74 

31 ± 109 

-63 ± 270 

782 ± 97 

584 ± 146 

1302 ± 386 

Si XHI 

42 ±64 

88 ± 191 

-60 ±21 

488 ± 69 

704 ± 361 

506 ± 79 

S XV 

99 ± 357 

1168 ± 1203 

-27 ± 633 

540 ± 206 

966 ± 1256 

69 ±1254 


1200 


1000 


' (/) 

E 800 

X 

X 600 


400 

500 1000 1500 2000 2500 

Excitation Energy (eV) 

Fig. 7.— Half-widths of the H-like Lya lines (km s“^) and the He-like resonance lines versus excitation 
energy (eV) of the upper level of the transition. The full and dashed lines represent the best linear ht to the 
HWHMs from the H-like lines, and the He-like lines (excluding the O VII width), respectively. 
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Figurej^shows the dependence of the half width 
at half maximum of the Gaussian fit versus the 
excitation energy of the upper level of the tran¬ 
sition. The linear correlation coefficient for the 
H-like half-widths is —0.89, indicating a strong 
anti-correlation between line half-width and ex¬ 
citation energy. For the He-like lines, the linear 
correlation coefficient is —0.81, also indicating a 
strong anti-correlation of line half-widths and ex¬ 
citation energy. Thus the line widths are anti¬ 
correlated with the upper energy level, in that 
the line width decreases with excitation energy. 
This anti-correlation shows that the more highly 
excited lines form at lower velocities, and thus 
closer to the stellar surface of the primary, indi¬ 
cating that the higher-temperature X-ray emission 
emerges from deeper regions in the wind than the 
cooler emission. 


In Figure the O VII line width seems lower 
compared to the trend defined by the more highly 
excited ions. Excluding the O VII line, a linear fit 
to the remaining He-like lines yields a linear cor¬ 
relation coefficient of —0.87, indicating a stronger 
anti-correlation, and also results in a steeper linear 
slope. This linear fit predicts that the O VII line 
should have a half-width of 918 eV, a factor of 1.2 
larger than observed. We caution that, unlike the 
other lines, the O VII line was only observed in one 
grating order since ACIS-S chip 5 was switched off 
during these observations. 


As a crude approximation, if we assume that 
the X-ray emitting material resides in a thin 
spherical shell at radius r around 6 Ori Aal, 
then the line prohle will extend from —V{r) to 
+V{r)^l - {RAai/rY, where RAai is the radius 
of 5 Ori Aal, and V{r) = Eoo.Aai(l - RAai/r)^, 


the standard velocity law for radiatively driven 
winds. The inverse correlation of the line widths 
with excitation energy suggests that the hotter 
X-ray emitting gas is formed over a smaller vol¬ 
ume in the wind acceleration zone closer to the 
star, where wind radial velocity differentials are 
larger and where higher temperature shocks can 
be generated; cooler ions can be maintained far¬ 
ther out in the wind where the acceleration (and 
thus the velocity differential) is smaller. A similar 
conclusion was reached by Herve et al. (2013) in 
their analysis of C Puppis. 


6.2.2. Effects of X-ray Line Opacity 


The possibility that strong resonance line pho¬ 
tons might be scattered out of the line of sight 
has signihcant implications on our physical un¬ 
derstanding of the X-ray emission from hot stars, 
especially in the interpretation of mass-loss rates 
derived from X-ray line profiles and abundances 
derived from X-ray line ratios. Resonance scatter¬ 
ing may be important for lines with high oscilla¬ 
tor strengths and could, in principle, change the 
line shape or intensity ratios, though recent anal¬ 
ysis by Bernitt et al. (2012) suggested that our 
poor knowledge of the underlying atomic physics 
may play the dominant role in accounting for 
discrep ancies in line intensities. Miller et al. 

(2002) focussed on the Fe XVII lines at 15.014 A 
and 15.261 A, which have oscillator strengths of 
2.49 and 0.64, respectively. Resonance scatter¬ 
ing might significantly affect the 15.014 A emis¬ 
sion line, which is one of the strongest lines in the 
5 Ori A X-ray spectrum, while scattering should 


be unimportant for the weak 15.261 A line. Miller 


et al. (2002) found that the observed ratio of these 


two lines, as derived from their Chandra grating 
spectrum, was /i5.oi/7i5.26 = 2.4 ± 1.3, nominally 
(though not signihcantly) below the optically thin 


limit /i5.oi/7i 5.26 = 3.5 derived from the Smith 
fc Brickhousej (2000) version of the Astrophysical 


Plasma Emission Code (APEC). 


We re-examined this issue for these two Fe XVII 
lines using our deeper spectrum and a slightly dif¬ 
ferent technique. We isolated the Fe XVII line 
region in the combined spectrum and fit this re¬ 
stricted region with an APEC-derived model, with 
abundances fixed at solar, including line broaden¬ 
ing. We first fit the Fe XVII line at 15.261 A, 
ignoring the region around the stronger 15.014 A 
line. We then included the 15.014 A line region 
and compared the predicted strength of the model 
15.014 A line to the observed line. This tech¬ 
nique, in which we use a full thermal model to 
fit the spectra rather than a simple comparison 
of line intensities, has the benefit that line blends 
in the region will be more properly taken into ac¬ 
count. We found that the model based on the best 
fit to the 15.261 A line greatly overpredicted the 
strength of the 15.014 A line, and can be ruled 
out at high conhdence {Xu = 3.57, restricted to 
the 14.90-15.14 A region; excluding this region, 
= 0.72). This may be an indication of the effect 
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of resonance scattering on the 15.014 A Fe XVII 
line. Since it appears that the 15.014 A line is a 
bit narrower than the 15.261 A line, we also re-did 
the fit, allowing the width of the 15.014 A line to 
differ from that of the 15.261 A line. We then re-fit 
only the 15.014 A line, allowing the line broaden¬ 
ing to vary and also allowing the normalization 
to vary. Figure shows the resulting fit. The 
best-fit HWHMs for the 15.014 A and 15.261 A 
lines are 1275l'l268 1496lit°3 km s“^, 

respectively, while the model normalizations are 
0.0024^°;°°?^ and 0.0030lJ^;!|J^i for the 15.014 A 
and 15.261 A lines, respectively. This analysis 
also shows the 15.014 A line is signihcantly weaker 
than expected compared to the 15.261 A line. This 
again may indicate that resonance scattering plays 
a role in determining the line profile shape and 
line strength, at least for the Fe XVII line, though 
uncertainties in the atomic models and in our def¬ 
inition of the temperature distribution for 6 Ori A 
may play a significant role in altering the intensity 
ratios for these lines. 

To further investigate the importance of reso¬ 
nance scattering, we also considered the Ne X lines 
at 10.239 A and at 12.132 A, which have oscillator 
strengths of 0.052 and 0.28, respectively. These 
lines complement the Fe XVII analysis since for 
Ne X the stronger line appears at longer wave¬ 
length; this means that any effects of differen¬ 
tial absorption that might affect the Fe XVII line 
analysis would have the opposite effect on the 
Ne X lines. We again fit the Ne X 10.239 A line 
with a single temperature APEC model, but hxed 
the temperature to the temperature of maximum 
emissivity of the Ne X lines, i.e. T = 6.3 x 10® K. 
We then compared the model that best fits the 
Ne X 10.239 A line to the Ne X 12.132 A line. 
Note that the Ne X 12.132 A line is blended with 
the Fe XXI line at 12.285 A (which has a tem¬ 
perature of maximum emissivity of 12.6 x 10® K, 
about twice that of the Ne X line), so we re¬ 
stricted the Ne X 12.132 A fitting region to the 
interval 12.0-12.22 A. We again find that the 
model, which provides a good fit to the weaker 
line {xl = 0.79), overpredicts the strength of the 
stronger line {xt = 8.63), again a possible indi¬ 
cation that resonance scattering is important in 
determining the flux of the strong line. 


7. The Influence of Colliding Winds on the 
Embedded X-ray Emission 

Colliding winds can have important observ¬ 
able effects in our analysis of the X-ray emis¬ 
sion from S Ori Aa in two ways. The collision 
of the primary wind with the surface or wind 
of the secondary could produce hot shocked gas 
which might contaminate the X-ray emission from 
the embedded wind shocks in the primary’s un¬ 
perturbed wind. In addition, the colliding wind 
“bow shock” around the weaker-wind secondary 
produces a low-density cavity in the primary wind, 
and this cavity, dominated by the weak wind of 
S Ori Aa2, should show little emission from em¬ 
bedded wind shocks. Along the line between the 
stars, the stellar winds will collide at the point 
at which their ram pressures pv'j_ are equal (e.g. 
Stevens et al.||1992 |. Using the stellar, wind, and 
orbital parameters in Table Figure shows the 
ram pressures for Aal (solid) and Aa2 (dashed: 
apastron, dotted: periastron) assuming that the 
wind from each star follows the standard /? veloc¬ 
ity law, V{r) = noo(l — where U(r) is the 

wind radial velocity at a distance r from the star, 
R is the stellar radius, and we assume that /3 = 0.8 
or 1.0. The ram pressure of Aal’s wind is greater 
than that of Aa2 throughout the orbit, so the wind 
from Aal should directly impact Aa2’s surface, in 
this simple analysis. 

A more thorough treatment includes the effects 
of Aa2’s radiation on the wind of Aal (and vice 
versa). These effects include “radiative inhibition” 
( Stevens fc Pollock|1994 ) in which Aal’s wind ac¬ 
celeration along the line between the stars is re¬ 
duced by Aa2’s radiative force acting in opposition 
to the wind flow, and “sudden radiative braking” 


(Owocki & Gayley 1995 Gayley et al. 1997), where 


Aal’s strong wind, which would otherwise impact 
the surface of Aa2, is suddenly decelerated by 
Aa2’s radiation just above the surface of Aa2. To 
estimate the magnitude of these effects, we solve 
the ID equation of motion along the line of centers, 
accounting for both star’s radiative forces via the 
standard Castor, Abbott, and Klein (CAK) line 


forces ( 

Castor et al. 

1975 

) including the finite disk 

correction factor ( 

Mend & Abbott 1986| 

Paul- 


drach et al. 19861 and gravitational acceleration. 


We determine the CAK parameters Q and a (Gay¬ 


ley 1995) to yield the desired mass-loss rates and 
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Fit to the Fe XVI115.26 line, and comparison to the 15.01 line 


Fit to the NeX 10.24 line, and comparison to the 12.13 line 




Fig. 8.— Left: APEC-based modeling of the Fe XVII 15.014 A vs 15.261 A lines. We first fit the 15.261 A 
line by itself. The thick histogram compares that model to the observed spectrum in the 14.5-15.6A range. 
This shows that the model tgat fits the 15.26 A line overpredicts the strength of the 15.0lA line. The vertical 
lines from the continuum to the X-axis at 14.9 A and 15.14 A show the adopted wavelength range of the 
15.014 A FeXVII line. Right: APEC model fit to the Ne X 10.24 A compared to the NeX 12.134 A line. 
The model (shown by the thick histogram) that fits the weaker 10.24 A line overpredicts the strength of the 
stronger 12.134 A line. 



Fig. 9.— Left: Ram pressure of Aal (solid) and Aa2 at apastron (dashed) and periastron (dotted). The 
black lines show a (5=1 law, while the gray lines show a /3=0.8 law. The gray vertical lines represent the 
location of Aa2’s surface for these two phases. Right: ID solution to the equation of motion of the primary 
wind along the line between the stars (solid) at three different separations — apastron (top), semi-major axis 
(middle), and periastron (bottom). For comparison, the dashed curve shows a /3=0.8 law, and the dashed 
line shows terminal velocity. 
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Fig. 10.— Density (left) and temperature (right) structure in the orbital plane from the SPH simulation of 

5 Ori Aal (larger star) and 5 Ori Aa2 (smaller star). The arrow shows the orientation of the line of sight. 
The system is pictured at phase (p = 0.87. The collision of the wind from 5 Ori Aal against S Ori Aa2 
produces a low-density cavity in the wind of <5 Ori Aal, where the emission from embedded wind shocks is 
reduced. The collision also produces a layer of hot shocked gas at the boundary of the cavity which produces 
< 10% of the emission from the wind shocks embedded in the unperturbed wind from 6 Ori Aal. In the 
temperature plot on the right, the hot gas from embedded wind shocks in the winds from 6 Ori Aal and 

6 Ori Aa2 is ignored, to emphasize the hot gas along the wind collision boundary. 



terminal speeds for each star by using the standard 
reduction in mass-loss rate from the finite disk cor¬ 
rection factor, i.e., Mfd = Mcak/(1 + 

We numerically integrate the equation of motion 
to distances far from the star to yield the terminal 
velocity. Then we repeat the process including the 
radiation and gravity of both stars to determine 
the speed of each wind along the line between the 
stars. 

Figure shows the equation-of-motion solution 
for the primary wind. The initial velocity corre¬ 
sponds to a /3 = 0.8 law, but radiative inhibition 
causes the wind (solid) to accelerate less compared 
to the unmodified /3-law (dashed). In addition, the 
primary wind velocity does begin to decrease from 
radiative braking. However, Star Aa2’s surface is 
located at the end of each line, so that the pri¬ 
mary wind does not completely stop before it im¬ 
pacts the secondary surface. This indicates that 
the wind from star Aal should still impact the 
surface of Aa2, even when the influence of the ra¬ 
diation field of star Aa2 is taken into account. Fur¬ 
thermore, due to the strong radiation of Aal, the 
wind of Aa2 does not accelerate off the surface of 


the star toward Aal, further suggesting that Aal’s 
wind will directly impact Aa2’s surface. 

We used a 3D smoothed particle hydrodynam¬ 
ics (SPH) code developed by Benz (|1990) and Bate 


et al. (19951 to model the effects of the wind-wind 


collision on the extended system wind. Okazaki 


et al. (2008) was the first to apply this code to 


a colliding-wind system, and Russell (2013) and 


Madura et al. (2013) describe the current capa¬ 


bilities of the code, which we briefly state here. 
The stars are represented as two point masses, and 
throughout their orbit they inject SPH particles 
into the simulation volume to represent their stel¬ 
lar winds. The SPH particles are accelerated away 
from their respective stars according to a /3=1 law 
(absent from any influence from the companion’s 
radiation) by invoking a radiative force with a ra¬ 
dially varying opacity K(r), i.e. g^ad = K(r)F/c, 
where F is the stellar flux. We take effects of the 
occultation of one star’s radiation by the other star 
into account. Radiative inhibition is included in 
the code (within the context of the radially vary¬ 
ing opacity method), but radiative braking is not 
since it requires the full CAK solution for the wind 
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driving, which is not yet included in the SPH code. 
Radiative cooling is implemented via the Exact In¬ 
tegration Scheme ( Townsendj2009 1, and the abun¬ 


dances of both winds are assumed to be solar (As- 
plund et al.||2009 |. 

The importance of radiative cooling of the 
shocked material is determined by the p arame- 
ter X = (ii 2 Wg/M _7 (Stevens et al. 1992), where 
di 2 is the distance to the shock in 10^^ cm, vs is 
the preshock velocity in 10®cms“^, and M-^ is 
the mass-loss rate in 10“^MQyr“^. y > 1 in¬ 
dicates adiabatic expansion is more important, 
while X < 1 indicates that the shocked gas will 
cool radiatively. For the /3=1 law, x ranges from 
0-5 ^ X ^ 1-3 between periastron to apastron, so 
the shocked gas should cool through a combina¬ 
tion of adiabatic expansion and radiation. 


Figure shows the density and temperature 
structure of the interacting winds in the orbital 
plane using the parameters in Table The pri¬ 
mary wind impacts the secondary star as expected 
from the analytical treatment above, where it 
shocks with newly injected secondary SPH par¬ 
ticles. If this interaction leads to SPH particles, 
either belonging to Aal or Aa2, going within the 
boundary of the secondary star, these particles are 
accreted, i.e. removed from the simulation. The 
temperature plot of figure [l0| shows that this leads 
to hot, shocked gas around Aa2, but this must be 
deemed approximate since the code does not force 
the Aal particles to accrete at the sound speed, 
which would increase the shock temperature, nor 
does it include any reflection of Aal’s radiation 
off of the surface of Aa2, which would decrease 
the shock temperature. The half-opening angle is 
^ 30°, so ^ 8% of the solid angle of Aal’s wind is 
evacuated by Aa2 and its wind. 


To determine the X-ray flux from the wind- 
wind/wind-star collision, we solve the formal so¬ 
lution to radiative transfer along a grid of rays 
through the SPH simulation volume, for which we 


use the SPH visualization program Splash (Price 


20071 as our basis. The emissivity is from the APEC 


model ( 

Smith et al. 

2001 

) obtained from XSPEC 

( Arnaud|1996 

), the circumstellar material absorbs 


according to the windtabs model (Leutenegger 
et al.|2010 |, and the interstellar absorption is from 


TBabs (Wilms et al. 2000bI. The radiative transfer 
calculation is performed at 170 energies logarith¬ 
mically spaced from 0.2 to 10 keV (100 per dex). 


and generates surface brightness maps for each en¬ 
ergy. These are then summed to determine the 
model spectrum, and finally folded through X-ray 
telescope response functions to directly compare 
with observations. The overall contamination level 
of wind-wind/wind-star collision X-rays is < 10% 
of the Chandra zeroth-order ACIS-S observation, 
so the influence of emission from shocked gas along 
the wind-wind boundary is not very significant , 
though contamination may be larger in some re¬ 
gions of the spectrum, depending on the emission- 
measure temperature distribution of the colliding- 
wind X-rays compared to that of the X-rays arising 
from embedded wind shocks. We caution, how¬ 
ever, that the model X-ray flux is dependent on 
the boundary condition imposed at the surface of 
Aa2, and so imposing a condition where the in¬ 
coming wind from star Aal shocks more strongly 
(weakly) will increase (decrease) the amount of X- 
ray emission from the wind-star collision. 


8. Conclusions 

Delta Ori Aa is an X-ray bright, nearby, eclips¬ 
ing binary and so offers the potential to directly 
probe the X-ray emitting gas distribntion in the 
primary star’s wind as the secondary star revolves 
through the primary’s wind. Our Chandra pro¬ 
gram was designed to obtain high signal-to-noise 
and high spectral resolution spectrometry of this 
system throughout an entire orbit. In this paper, 
we have sought to characterize the overall spec¬ 
trum at its highest signal-to-noise ratio by combin¬ 
ing all of the Chandra spectra and examining tem¬ 
perature distributions and line parameters. Our 
main results are presented below. 


1. Our analysis of the Chandra image shows 
that the emission is mostly dominated by 
5 Ori Aa, with little detectable emission 
from S Ori Ab. 


2. The temperature distribution of the X-ray 
emitting gas can be characterized by three 
dominant temperatures, which agrees fairly 
well with the temperature distributions de¬ 
rived by the earlier analysis of|Miller et al. 


(2002) and Raassen & Pollock (2013). 


3. The strong lines are generally symmetric, 
and Gaussian profiles provide a reasonable 
representation of the profile shape, though 
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in most cases, and especially for the Ne X 
and Fe XVII there are significant deviations 
from Gaussian symmetry. 


4. The line widths determined by Gaussian 
modeling shows that half-widths are typi¬ 
cally 0.3—0.5x Foo, where Foo is the terminal 
velocity of the wind of i5 Ori Aal. These val¬ 
ues are generally larger than the line widths 
measured by Miller et al. (20021, though it is 
unclear whether this represents a real change 
in the line profile or if there is a calibration 
issue in the analysis of the earlier data set, 
which was obtained at an anomalously high 
focal plane temperature. 


5. We hnd a strong anti-correlation between 
the widths of the H-like and He-like tran¬ 
sitions and the excitation energy. This indi¬ 
cates that the lower-energy transitions occur 
in a region with larger velocities. Assuming 
a standard wind acceleration law, this cor¬ 
relation probably indicates that the lower- 
energy lines emerge from further out in the 
wind. 


6 . Analysis of strong and weak transitions of 
Fe XVII and Ne X indicates that resonance 
scattering may be important in determining 
the flux and/or shape of the stronger line. 
This agrees with the analysis of the Fe XVII 
line by Miller et al. (20021 but at higher sig¬ 
nificance. We caution that some of these 
differences in the observed to predicted line 
ratios may be influenced by an inaccurate 
temperature distribution and/or uncertain¬ 
ties in the atomic physics. It is also inter¬ 
esting to note that these two lines also have 
the most non-Gaussian profiles, as shown in 
Figure perhaps indicative that some line 
photons have been scattered out of the line 
of sight. 


The spectrum combined from the four indi¬ 
vidual G/iandra-HETGS observations represents a 
very high signal-to-noise view of the emission from 
<5 Ori Aa. However, these observations were ob¬ 
tained at a variety of orbital phases, so that the 
combined spectrum is a phase-averaged view of 
the overall X-ray emission from S Ori Aa. In a 
companion paper (Nichols et al.|2015 1 we look for 


the effects of phase- and time-dependent changes 
in the continuum and line spectrum. 
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